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In the standard model, stabilization of a classically unstable cosmic string may occur through the 
quantum fluctuations of a heavy fermion doublet. We review numerical results from a semiclassical 
expansion in a reduced version of the standard model. In this expansion the leading quantum 
corrections emerge at one loop level for many internal degrees of freedom. The resulting vacuum 
polarization energy and the binding energy of occupied fermion energy levels are of the same order, 
and must therefore be treated on equal footing. Populating these bound states lowers the total 
energy compared to the same number of free fermions. Charged strings are already stabilized for a 
fermion mass only somewhat larger than the top quark mass. Though obtained in a reduced version 
these results suggest that neither extraordinarily large fermion masses nor unrealistic couplings are 
required to bind a cosmic string in the standard model. Furthermore we also review results for a 
quantum stabilization mechanism that prevents closed Nielsen-Olesen type strings from collapsing. 

PACS numbers: 11.27.+d, 03.65.Ge, 14.65.Jk 


I. INTRODUCTION 

In this brief review article we summarize some recent results on how quantum effects can stabilize classically unstable 
string-like configurations in the electroweak standard model. Such configurations are the particle physics analogs of 
vortices or magnetic flux tubes in condensed matter physics. They are usually called cosmic strings to distinguish 
them from the fundamental variables in string theory, and also to indicate that they typically stretch over large 
length scales. For electroweak strings or Z-stringslIHS], the Z-component of the electroweak gauge boson acquires 
the structure of an Abelian Nielsen-Olesen type vortex [3]. 

Such strings may have emerged copiously at various epochs in the early universe, at interfaces between regions 
of different vacuum expectation values of the Higgs field(s) before electroweak symmetry breaking was reached. If 
these strings are stable, they should have survived and we should be able to observe them today. While their 
direct gravitational effects are negligible, Z-strings can still be relevant for cosmology at a sub-dominant level [SJ E]. 
Their most interesting consequences originate, however, from their coupling to the standard model fields. Z-strings 
provide a source for primordial magnetic fields [3] and they also offer a scenario for baryogenesis with a second 
order phase transition [71 [5] ■ In contrast, a strong first order transition, as required by the usual bubble nucleation 
scenario, is almost certainly ruled out in the electroweak standard model|9| without non-standard additions such as 
supersymmetry or higher-dimensional ODerators|10j. For reviews on cosmic strings and cosmological implications of 
networks of strings see e.g. Refs. [IIHS]. 

There are numerical simulations showing that distributions or networks of intersecting strings are stable and ap¬ 
proach a so-called scaling limit, cf. Ref. [13] and references therein. However, the standard model has neither a 
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classical nor a topological mechanism to stabilize an isolated string to ensure that we would still observe string net¬ 
works at the present epoch. Thus these interesting configurations are only viable if they are bound dynamically by 
their interaction with the quantum fluctuations of the fields. The most important contributions are expected to come 
from (heavy) fermions, since their quantum energy dominates in the limit Nc —> oo, where Nc is the number of any 
other internal degree of freedom such as QCD color, for instance. The Dirac spectrum in typical string backgrounds 
is deformed to contain either an exact or near zero mode, so that fermions can substantially lower their energy by 
binding to the string. This binding effect can overcome the classical energy required to form the string background. 
However, the remaining spectrum of modes is also deformed and for consistency its contribution (the zero-point or 
vacuum polarization energy) must be taken into account as well. Heavier fermions are expected to provide more 
binding since the energy gain per fermion charge is higher and the Yukawa coupling is stronger. A similar scenario is 
also suggested by decoupling arguments [TS]. 

Strings coupled to quantum fields have been investigated previously. Naculich|lt)j has shown that fermion fluc¬ 
tuations destabilize the string in the limit of weak coupling. The quantum properties of Z-strings have also been 
connected to non-perturbative anomalies [T7]. The emergence or absence of exact neutrino zero modes in a Z-string 
background and the possible consequences for the string topology were investigated in Ref. [TB]. A first attempt at 
a full calculation of the fermionic quantum corrections to the Z-string energy was carried out in Ref. |19j . However, 
those authors were only able to compare the energies of two string conhgurations because of limitations arising from 
the non-trivial behavior at spatial infinity. We will present a solution to this serious obstacle. The fermionic vacuum 
polarization energy of the simpler Abelian Nielsen-Olesen vortex has been estimated in Ref. [20] with a non-standard 
regularization that only subtracts the divergences within the heat-kernel expansion. Quantum energies of bosonic 
fluctuations in string backgrounds were also calculated in Ref. m- Finally, the dynamical helds coupled to the string 
can also result in (Abelian or non-Abelian) currents running along the core of the string. The time evolution of such 
structured strings was studied in Ref. [52], where the current was induced by the coupling to an extra scalar field. 

This review is based on a number of publications [231j26] studying the contribution of fermion quantum corrections to 
the vacuum polarization energy of a straight and infinitely long cosmic string. The technical details of the underlying 
calculation are thoroughly described in the appendices of Ref. |2S] • 

Cosmic strings normally carry magnetic flux and must either extend to spatial infinity or be closed. In the latter 
scenario, the straight strings described above would be replaced by quantum stabilized toroidal strings, or torons. 
A simple application of Derrick’s theorem m reveals that torons are classically unstable, but this conclusion may 
again be reversed by quantum corrections. In section VI a simple application[28] of Heisenberg’s uncertainty principle 
suggests that quantum fluctuations indeed prevent toroidal strings from collapsing. 

This review is organized as follows: In sectionjlTjwe introduce the model that we use to investigate the interaction of 
a cosmic string with a fermion doublet. In section |TTl| we explain the essential ingredients of the spectral methods that 
allow us to compute the vacuum polarization energy for string-like conhgurations that are translationally invariant 
with respect to one or more coordinates. We present numerical results for the vacuum energy of a string background 
in section jlV] and discuss the binding of charged strings in sections jVj In section jV^ we briefly report some novel 
results based on enlarging the limited subset of parameters used in the previous numerical analysis. In section VH we 
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review a quantum mechanical mechanism to stabilize closed Nielsen-Olesen type strings finally we conclude in section 

ivnii 


II. MODEL 


For the current investigation the fermion doublet will be considered degenerate so that the introduction of a matrix 
notation for the Higgs field is appropriate. Two angles and ^2 that parameterize the isospin orientation|29| serve 
as variational parameters of the string configuration. The isospin group SU{2) can, in general, be described by three 
Euler angles; in the considered configurations the third angle picks up the winding of the string in azimuthal direction. 
For notational simplicity, we introduce the abbreviations Si = sin(.fi) and = cos(^i). Then the string configuration 
reads 


= vlnip) 


/ S1S2 e-“‘^ 
\-ici + S 1 C 2 


—ici — S 1 C 2 
S 1 S 2 


( 1 ) 


for the Higgs field and 


W = W“ — = n Si S 2 & 

2 gp 




' S1S2 - (fci + S1C2) e™”^^ 

(ici - S1C2) -S1S2 , 


( 2 ) 


for the gauge boson (in temporal gauge).^ The variables p and ip are polar coordinates in the plane perpendicular to 
the string, while the Higgs vacuum expectation value v and the gauge coupling constant g are model parameters. The 
string configuration involves profile functions fn and fc which are the analogs of the Nielsen-Olesen vortex profiles 
with boundary conditions 


/g , /ff -t 0 

for p —0 

(3) 

/g , /ft -t 1 

for p —>■ 00 . 

(4) 


The integer n is the winding number of the string, for which we will typically take n = 1 in the numerical calculations. 

This simplified version of the standard model has a vanishing Weinberg angle so that the t/(l) hypercharge decouples 
and the gauge symmetry is SU{2)l. Then the (classical) boson part of the Lagrangian reads 

Cbos = -^tr (G'^'G^,) + itr (D^$)^ D^<i> - ^tr ($^$ - v^)\ (5) 

with the covariant derivative — igW^ and the SU{2) field strength tensor 

G^, = d^W,- - ig [W^, W,] . (6) 

The classical boson masses are determined from g and v and the Higgs self-coupling A as mw = gvjV^ and mn = 
2v y/\ for the gauge and Higgs bosons, respectively^. The interaction with the fermion doublet is described by the 


^ The gauge field is a vector both in coordinate and isospin space, and the generators of the isospin algebra are (t“/ 2 with the Pauli-matrices 
^ The three gauge bosons are degenerate due to the vanishing Weinberg angle. 
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Lagrangian 

Aer = {PlP+ PrP) ^ ^ {^Pr + $^Pl) 4- (7) 

with the left/right-handed projector Pr, i = (il ± 75 ) /2. Upon spontaneous symmetry breaking the Yukawa coupling 
/ induces a fermion mass m = vf. Assuming a heavy fermion doublet with the mass of the top quark, the standard 
model suggests the parameters 


5 = 0.72, u = 177GeV, TOH = 140GeV, / = 0.99. 


( 8 ) 


In the numerical search for a stable string, other model parameters are considered as well, cf. section 
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The classical energy per unit length of the string is determined by £bosj 
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where the dimensionless radial integration variable is related to the physical radius by Pphys = p/r^i ^md we have 
introduced the mass ratio fiH = mYilm. 

The central object of the present investigation is the fermion contribution to the energy. It is obtained from the 
solutions to the Dirac equation in the two-dimensional plane perpendicular to the string (in a basis with 70 diagonal) 


i74'„ = 


with 




VO" ■ v? 


a-(P 

0 ) 


d^ + m 



+ 77int , 


( 10 ) 


where the single particle Hamiltonian H is extracted from £fer, eq. Q. Details of the interaction part, Hint, will be 
discussed later, in particular concerning the subtleties at spatial infinity. The spectrum of the Dirac operator consists 
of bound state solutions with discrete eigenvalues a;„ = and continuous scattering solutions whose eigenvalues w 
are labeled by momentum k, i.e. oj = y/k"^ + m?. Note also that the Dirac Hamiltonian for the string background 
anti-commutes with = pp so that the spectrum is charge conjugation invariant and it suffices to consider the 
non-negative eigenvalues. 


III. FORMALISM 


The vacuum polarization energy is the renormalized sum of the changes of the zero-point energies of fermions^ in 
the background of a static configuration 


F --- 

-^vac — 2 




h 


= dk uJk PrePk). 


( 11 ) 


Here a;„ are the energy eigenvalues in the presence of the string as obtained from the Dirac equation (10), and p 


,(o) 


are their free counterparts. In the second part of equation (11 1 the changes of the single particle energies are split up 


into the contribution from distinct bound states (e^) and an integral over the energies of the continuous scattering 


® For bosons the overall sign needs to be reversed. 
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states weighted by the (renormalized) change of the density of scattering states, A p^enik). This change is induced by 
the background helds and its determination will be outlined in the next subsection. The scattering states obey the 


standard dispersion relation ojk = \/k'^ + In eq. (11), the factor h has been made explicit to stress that Svac is 
a quantum effect. Unless stated otherwise, natural units (h = 1 and c = 1) will be adopted in the remainder of this 
paper. 

The computation of the fermion contribution to the vacuum polarization energy of a cosmic string proceeds in three 
stages. First, spectral methods are employed to express this energy in form of scattering data|30j while maintaining 
standard renormalization conditions. Second, this approach is extended to accommodate configurations that are 
translationally invariant in one or more spatial directions, using the so-called interface formalism |31j. Third, the 
non-trivial structure of the cosmic string configuration at spatial infinity must be addressed. More precisely, since 
the fields approach a pure gauge rather than zero at spatial infinity, the immediate application of spectral methods 
is impossible. Instead, a particular local gauge transformation is performed, which facilitates the formulation of a 
well-behaved scattering problem [24]. 


A. Spectral Methods 


The basic idea of the spectral approach is to express Ap^g^ik) as the momentum derivative of the phase shifts of 
the scattering states [35]. Then 


K, 






( 12 ) 


In this formulation, it is implicitly assumed that the system has enough symmetry to allow for a partial wave 
decomposition with a degeneracy factor Di for the partial wave of (generalized) angular momentum £. The major 


concern in eq. (12) is renormalization. As it stands, the vacuum polarization energy eq. (12) in three spatial dimensions 


is quadratically divergent at large momenta. In this regime, however, the Born series adequately represents the phase 
shift. It is therefore advantageous to introduce 




( 13 ) 


where the subscript on the phase shift indicates the subtraction of the N leading terms of the Born series to the phase 
shift, = Se — — ... — Choosing N sufficiently large renders finite. 

Technically the Born series is an expansion in powers of The contribution of any such power to the vacuum 

polarization energy can be associated with a Feynman diagram 




27r dk 



(14) 


where the loop corresponds to the quantum fluctuation, and the external lines represent insertions of the background 


field, i.e. the interaction of the cosmic string with the fermions. Therefore, the Born subtraction in equation (13) 


corresponds to the sum E^ of all Feynman diagrams with up to N insertions of the background interaction. These 
diagrams can be computed with standard techniques; in particular dimensional regularization can be implemented 
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to handle the ultra-violet divergences. It must be emphasized at this point that the use of the Born series does not 
imply any approximation for the phase shift; it is always computed exactly (eventually numerically). The Born series 
merely serves to identify the divergences at large momenta in the spectral approach. 

The main advantage of transferring of the divergences into Feynman diagrams is that they can now be straight¬ 
forwardly combined with the counterterm contribution to the energy, i?cT- The latter is found by substituting the 
background configuration into the counterterm Lagrangian in just the same way as the classical energy, eq. § is 
obtained from C^os ■ In any (multiplicatively) renormalizable theory, the counterterm Lagrangian has the same struc¬ 
ture as the classical one and a suitable choice of its coupling constants cancels all ultra-violet divergences completely. 
Then + ifcT is free of any divergences. The finite pieces of the counterterm coupling constants will be uniquely 
determined from appropriate conditions describing properties of the particles that are associated with the fields. This 
procedure will be briefly discussed in section. |IV C| In total, the sum 

i^vac = E^^^ + + Ect (15) 


gives an unambiguous result for the vacuum polarization energy once the renormalization conditions are fixed. A first 
principle derivation of this result, based on a quantum field theoretic formulation of the energy momentum tensor and 
the analytic properties of the Greens function, is presented in Ref. [33] • We stress that eq. (15) is a sum of separately 
finite contributions, which does not involve any large number {e.g. a cut-off). 


B. Interface Formalism 


Equation (15) is perfectly suited for a computation of the vacuum polarization energy for static configurations that 


allow for a full decomposition into partial waves. However, the cosmic string does not exhibit a three-dimensional 
rotational invariance; rather it is translationally invariant along the z-axis. In this scenario the wave-function of the 
quantum fluctuation factorizes into 






(16) 


Here is the reduced wave-function for a two-dimensional scattering problem in the plane perpendicular to the 

string with the total dispersion relation uj = + nrfi = -I- With the replacements —>■ e| and 

+ uj'^ in equation (^ISl, an integration over p (with measure dp/2TT) yields the vacuum polarization energy 


Wfc 


per unit length of the string. However, this procedure runs into an immediate problem: The scattering data do not 
depend on the momentum p and hence the p-integral will not be finite for any N. A careful analysis must treat the 
p-integral in dimensional regularization]^ 

2(47r)^ t [i Jo 


eW 


rlk ti+1 d 


N 


(17) 


where d is the analytic dimension of the subspace in which the configuration is translationally invariant. The divergence 
now manifests itself via the singularity of the T-function coefficient as d —> 1. Due to sum rules for scattering data[34j . 
which represent generalizations of Levinson’s theorem |35j . the expression in curly brackets in eq. 0 vanishes as 
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d —>■ 1 which cancels the F-function pole. Hence the limit d —>■ 1 can indeed be taken[5T]. The result is 


1 

47r 


f dk 

ojI In 


[ Jo TT 


kkJ J 


dk 


E 


In^ — 


(18) 


where is expressed through the non-negative part of the spectrum thanks to charge conjugation invariance in the 
string background. Here Hr is an arbitrary renormalization scale that has no effect on by exactly the same sum 


rules. The expression (181 for E^^'> replaces the analog in equation (^15|. When comparing the two expressions, it is 


apparent that the function multiplying the (Born subtracted) phase shift is of higher power in k than its counterpart 
before integrating over the momentum p. Hence N must be increased when further dimensions are added in which 
the configuration is translationally invariant. This, of course, merely reflects the fact that ultra-violet divergences 
become more severe in higher dimensions. 

The expressions obtained so far for the vacuum polarization energy rely strongly on the analytic properties of the 
scattering data. These properties also allow for numerical calculations using purely imaginary momenta. This has (at 
least) two advantages: (i) the oscillating phase shifts turn into exponentially decaying (logarithms of) dost functions, 
and (ii) the momentum integral and the sum over angular momenta may be exchanged. While (i) drastically improves 
numerical stability, (ii) significantly simplifies the treatment of the logarithmic divergences that emerge at the third 
and fourth order of the Born and Feynman expansions. These technical details are discussed at length in Ref. [26]; 
we address them briefly in section EnDj 


C. Choice of Gauge 


The results presented thus far establish a formalism for computing the renormalized vacuum polarization energy 
of background fields in the string geometry. However, there is stillthe problem that the string does not induce a 


well-behaved Born series when expanding in powers of iFint from eq. (10). Even though the full Hamiltonian is gauge 
invariant, iJint is not and does not vanish at spatial infinity. This obstacle appears because the Dirac Hamiltonian 
that is obtained by straightforward substitution of the string background, eqs. ([^ and ([^, does not turn into the 
free Dirac Hamiltonian at p ^ oo. Instead it becomes H —>■ W{(p)H{,;eeU{ip). This local gauge transformation is a 
consequence of the string winding and acts only on the left-handed fermions 


( S2COs{nip) \ 
—S 2 sin(n(/7) 
C2 J 


(19) 


Unfortunately, the obvious gauge transformation H —>■ U{ip)HW {(p) does not solve the problem completely: Although 
it would generate vanishing interactions at infinity, it will also induce a 1/p^ potential at the core of the string, p —)■ 0. 
This might still yield well-defined phase shifts, though they would likely be difficult to compute numerically. In any 
event, the conditions underlying the analyticity of the scattering data are certainly violated by this singular behavior. 
As argued at the end of the previous section, analyticity is central for numerical feasibility of our approach. As a 
solution, we can define a radially extended gauge transformation 


U{p,(p) = PlBxp {in ■ f ^{p)) + Pn . 


( 20 ) 










This transformation fixes the gauge and leads to the interaction term in equation (10) 


TJint = m 


(///cos(A) - 1) 


1 0 
0 -1 


+ ifH sin(A) 


0 1 
-1 0 


1 ' P S ■ p 

2dp\^-p -a ■ p 


n ■ T 


+ ■ 


ns2 / —a ■ (p u ■ (p 


2p 


a ■ ip —a ■ p 


fc sin(A) /g(A) + (/g - 1) sin(^(p)) /g(-C) 


The new gauge function ^{p) appears via the difference A(p) = — ^(p) and the isospin matrices are 


C 2 


S 2 e 




S2e-“‘^ -C2 


and Ig{x) = 


—S2sin(a;) 


[c 2 sin(a;) — icos( 2 :)] 


( 21 ) 


( 22 ) 


[c 2 sin(x) + icos(x)] e S 2 sin(a;) / 

Imposing the boundary conditions ^(0) =0 and C(oo) = for the new gauge function ^(p) defines a well-behaved 
scattering problem. The specific form of ^(p) is irrelevant (apart from its boundary conditions) and must not have 
any influence on the quantum energy, since it merely parameterizes a gauge transformation. 

Note that the gauge transformation is single-valued at spatial infinity, U{oo,p) = [/(oo, tp + 27r). In this respect it 
differs from the analogous problem of fractional fluxes in QED. In that case a similar choice of gauge would not be a 
remedy; rather the calculation of the vacuum polarization energy requires the introduction of a return flux to arrive 
at a well-behaved scattering problem [36]. The return-flux approach can also be used for the present calculation, but 
it is much more laborious numerically [23l|24]. 


D. Fake Boson Field 


The idea of utilizing a fake boson field to simplify the treatment of higher order divergences was first implemented 
in Ref. [38] . The third and fourth order contribution from the Born series produce logarithmic divergences but the 
corresponding Feynman diagrams are very cumbersome to evaluate. On the other hand, these logarithmic divergences 
are similar to the ones found in the second order vacuum polarization energy of a boson field fluctuating about a scalar 
potential. Matching its strength appropriately and recalling that for imaginary momenta the momentum integral and 
the angular momentum sum may safely be exchanged allows to replace 





d 

dt 


De (vi{t) - - Y 

. e e . 


(23) 


under the integral in eq. (18). The quantity with this replacement will be called Es- (Overlined quantities refer 


to the bosonic scattering data). The replacement eq. (23) must, of course, be accompanied by the boson Feynman 
diagram E^ so that the total vacuum polarization energy becomes 


i^vac = Es + + E\ 


(24) 


where the superscript indicates the inclusion of the counterterm contributions. Fach of the three terms on the right 


hand side of equation (24) is separately ultra-violet finite. The advantage of eq. (24) and the replacement eq. (23) is 
now obvious: Instead of fermionic contributions up to order TV = 4, we only need to compute second order fermionic 
and bosonic Feynman diagrams and terms in the corresponding Born series. 
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IV. NUMERICAL RESULTS FOR THE VACUUM POLARIZATION ENERGY 


Before discussing numerical results it is worthwhile to return to the issue of the particular gauge transformation, 


eq. (20). The resulting invariance provides an excellent reliability test for the numerical procedure: Modifying the 


shape of ^(p) while keeping its boundary conditions fixed alters Es and in eq. (24). However, these individual 


changes are bound to compensate each other [23]. We discuss this test on the numerical results in section B below. 

For ^2 f the Hamiltonian acquires imaginary parts, though it still remains hermitian. This complicates the 
computation of the phase shifts because in-coming and out-going spherical waves obey different differential equations. 
Therefore most of the numerical results are limited to the case ^2 = f as in ref. [25] . We will present some novel results 
without that restriction at the end of this section. 

As a convention, all numerical results and dimensionful parameters from here on will be measured in appropriate 
units of the (perturbative) fermion mass m. 


A. Variational Ansatze 


Despite the simplification in eq. (24), the numerical computation is still expensive. The scattering data are extracted 
from a multi-channel problem and for the final result to be reliable several hundred partial wave channels must be 
included. This numerical effort restricts the number of variational parameters that can be used to characterize the 
profile functions. In addition to ^ 1 , we introduce three scale parameters wh, ww and via the ansatze 

[l-e-(^)^ 


/ff(p) = l-e , /G(d) = l-e 


?(p) = Cl 


(25) 


The scale parameterizes the shape of the gauge profile which should not be observable c.f. the remark at the 
beginning of this chapter. The other properties of the profiles are chosen to keep E^i regular. 

We have also considered an exponential parameterization for the gauge field 

/g(p) = 1 - + —) exp (" —(26) 

V wgJ \ wgJ 

which yields a slightly better agreement with the original Nielsen-Olesen profiles which minimize E^i for = tt/2. 
No significant difference in ifyac was found between these ansatze. 


B. Gauge Invariance 

We check gauge invariance by varying the width of the gauge profile, C(p)- A typical result is shown in table As 
expected, the individual contributions to i?vac depend strongly on w^. However, these changes essentially compensate 
each other. Numerically the most cumbersome part of the calculation is Es- From various numerical tests (changing 
the extrapolation scheme for partial wave sum, the momentum integration grid, etc.) its numerical accuracy is 
estimated to be at the 1% level. Within that precision range, E^ac is independent of w^, which verifies both gauge 
invariance and the equivalence of terms in Born and Feynman series. This precise equivalence of separately divergent 
quantities is fundamental to the spectral approach. 
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TTren. 

-^FD 

Es 


Eyac 

0.5 

-0.2515 

0.3489 

0.0046 

0.1020 

1.0 

-0.0655 

0.1606 

0.0032 

0.0983 

2.0 

-0.0358 

0.1294 

0.0038 

0.0974 

3.0 

-0.0320 

0.1235 

0.0056 

0.0971 

4.0 

-0.0302 

0.1193 

0.0080 

0.0971 


TABLE I: Numerical results for the various contributions (241 to the fermion vacuum polarization energy in the minimal 
subtraction scheme. 


C. On—Shell Renormalization 

The results reviewed in the previous subsection were obtained in the MS renormalization scheme, which essentially 
omits the non-divergent parts of the Feynman diagrams. In this scheme, the dependence of F^vac on the model pa¬ 
rameters factorizes which simplifies the computation because this dependence can easily be traced from Ed- As a 
consequence, many parameter settings can be studied by rescaling instead of recomputing i?vac- Any other scheme 
merely differs by manifestly gauge invariant (finite) counterterms. Since physically meaningful results require renor¬ 
malization conditions that correspond to a particle interpretation, an additional (mild) parameter dependence in Afd 
is induced. 

To be specihc, we consider the so-called on-shell scheme, in which the coefficients of the four allowed counterterms 
(three terms from eq. ^ supplemented by tr $]) are determined such that 

• the tadpole graph vanishes implying that the Higgs vev v remains unchanged, 

• the Higgs mass remains unchanged 

• the normalization of one-particle states of the Higgs boson remains unchanged 

• the normalization of one-particle states of the vector meson remains unchanged 

in the presence of fermionic quantum corrections. Note that the vector meson mass M\y is not hxed by these conditions 
and thus will be a prediction that includes quantum corrections. This suggests to tune the gauge coupling to reproduce 

I 

an appropriate physical value for Mw= Mz ~ 90GeV. Typical results for the vacuum polarization energy per unit 
length of the string are shown in figure [l] as functions of the variational parameters. Except for narrow string 
configurations dominated by the Higgs field, the vacuum polarization energy turns out to be positive. Therefore, 
fermionic vacuum fluctuations alone do not provide any substantial binding and no stable uncharged string is found 
for the physically motivated parameters, eq. ([^, for which Ed dominates the total energy. Since Ed decreases quickly 
with increasing Yukawa coupling /, some stability is indeed seen for large / and narrow strings but in this regime 
the restriction to one fermion loop in the vacuum polarization energy is unreliable because of the unphysical Landau 
ghost contribution [saiini- 
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V. CHARGED STRINGS 


Cosmic strings induce many fermionic bound state levels for the two-dimensional scattering problem, whose energies 
are denoted by Cj. For = ^2 = 7’'/2 there even exists an exact zero mode[TB]. In the three-dimensional setting these 
bound states acquire a longitudinal momentum for the motion along the symmetry axis, i.e. their energies become 

/ T17T 

Ei{pn) = Y ef +pI with = —. (27) 


Here L is the length of the string. To leading order of the limit L —> 00 , the sum over the discrete longitudinal 
momentum turns into a continuum integral, ^ dp. To find the minimal bound state contribution at 

a fixed charge (per unit length), a chemical potential /i with max(|ej|) < /i < m is introduced, and all levels with 
Ei(p) < p. are populated. This procedure defines a Fermi momentum for each level, pf (p) = p^ — which enters 
the total charge per unit length of the string 

QW = — - ( 28 ) 

TT 

i 


This relation can be inverted to give p = p{Q) and thus pf = pf (Q). This quantity provides, for a given charge Q, 
the maximal momentum along the string direction for each fermion to get trapped in a two-dimensional bound state 
Ci perpendicular to the string. From this the binding energy (per unit length) for a prescribed charge 

rpf(Q) 


Ehind{Q) = - ^ 

TT ' ^ 


dp 


€-1 P‘^ — m 


(29) 


is computed relative to an equal number of free fermions that have energy m each. 

Figure shows the fermion contribution to the binding energy, i?vac + £'bind(Q)- In the left panel the graph 
terminates for a given configuration when all available bound state levels are occupied and the charge cannot be 
increased any further. For small charges, narrow strings are favorable while the binding energy of wider strings 
decreases more quickly as Q increases. Surprisingly, the envelope along which Flvac + F'bind is minimal (when varying 
the string parameters at fixed charge) forms a straight line. Extrapolating this line to Q = 0 indicates that the 
fermion vacuum polarization energy should (approximately) vanish. This extrapolation circumvents the Landau ghost 
problem mentioned in the previous section. 




FIG. 1: Fermion vacuum polarization energy. Left panel: on-shell renormalization scheme for configurations with ^2 = f. 
Right panel: some first results for arbitrary ^2 in the MS scheme. 
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^j=0.47Z 




FIG. 2: (Color online) Fermion contribution to the energy of a charged string (left panel). Total energy of the charged string 
(right panel) as function of the charge density. 


To finally decide on dynamical stability, the classical energy must be included as well. To this end several hundred 


configurations, that are characterized by the variational ansatze, eq. (25), were scanned. It is appropriate to label 
them by s = 1, 2 ,... and to compute their total binding energy as 


4:^(g) = +Nc te+ 


(30) 


for a given charge. If 


i^tot(Q) = min^ £'tot(<3) 


< 0 


(31) 


a stable configuration has been found. The right panel of figure shows iftot as a function of charge for various 
values of the Yukawa coupling constant. Since the mass m of a non-interacting fermion sets the scale in the numerical 
analysis, this refers to different values of the dimensionless ratio v/m. For / « 1.6 the classical and fermion energies 
essentially cancel each other and leave iftot roughly independent of charge"^. Bound objects are observed by further 
increasing the Yukawa coupling to about / « 1.7, which corresponds to a heavy fermion mass which is still less than 
twice the top quark mass. The minimizing configurations have « 0, i.e. they are dominated by the Higgs field. 
The numerical results reviewed here largely support the existence of a novel solution within the standard model of 
particle physics. These string configurations exist provided that a sufficient density of very heavy fermions is present 
which can be trapped to the string core to stabilize the configuration. Such conditions may have existed in the early 
universe, and the emerging cosmic strings should still be traceable today. Inverting this argument, not observing 
cosmic strings indicates that fermions coupling to the weak interactions in the usual way cannot have existed with 
sufficient density if their mass exceeds about twice the top-quark mass, unless they can decay to lighter fermions. 


4 


The strong dependence at small Q is artificial because very narrow strings have not been considered to avoid the Landau ghost 
inconsistency. 
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VI. VARIATIONAL PARAMETER ^2 


So far the numerical analysis of the charged string was limited to the case of ^2 = 7J'/2, for which the scattering 
and bound state problems only involve real valued Hamiltonians. First numerical results for the vacuum polarization 
energy at ^2 ^ l‘^ are shown in the right panel of figure They suggest that reducing ^2 may indeed further decrease 

the vacuum polarization energy. However, the contribution from explicitly occupying the bound fermion levels may 
have the opposite effect. The numerical determination of the bound state energies is more complicated when the 
Hamiltonian is not real. Rather than doubling the diagonalization problem by separating real and imaginary parts, 
it is easier to extract the bound state energies from the roots of the dost function at complex momentum. Since 
this dost function results from a multi-channel problem, its roots need not be simple which adds to the difficulty of 
extracting them numerically. Fortunately, Levinson’s theorem |34[ 135) yields the number of expected bound states 
from the phase shifts. The typical behavior of the bound state energies as a function of ^2 is shown in figure]^ Note 
that the zero mode for Ci = f and ^2 = f is recovered. Similarly to decreasing ^ 1 , any reduction of ^2 also reduces 
the binding of that level. 


.v.^.=/.a 4,=o.5n; 



..„=2.5. ...^=2.0, ^,=0.3n; 



FIG. 3: The variation of the zero modes as a function of ^2 for narrow (left panel) and wide (right panel) string profiles. The 
dots indicate calculated values while the lines have been added to gnide the eye. 

The data from the ^ 2 -variation can only lower the minimal total energy which further reduces the bound on / and 
thus the mass of the trapped fermion below which a charged cosmic string is stable. The conclusions made at the end 
of the previous section remain valid, but the bounds and statements become more stringent. 

VII. CLOSED NIELSEN-OLESEN STRINGS 

Cosmic strings usually carry magnetic flux and are thus particle analogs of Abrikosov vortices in type-H supercon¬ 
ductors. Since the flux is conserved, the vortex must either end on (magnetic) charges or be closed. The preferred 
topology of strings formed during the electroweak phase transition are thus closed rings which may be stabilized 
dynamically - either by spinning at high angular momenta m, by extra currents running along the ring [min], or 
by trapping charged particles as discussed earlier. Such stable vortex conhgurations are called vortons |43j . 

If none of these stabilizing effects is operative, the vorton will follow its (classical) instability and collapse by 
radiating away its energy. Again, this picture may be altered by quantum fluctuations: Due to the uncertainty 
principle, there is a quantum mechanical penalty for localizing in a very small volume. This quantum energy may 
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ultimately balance the classical energy of the vorton and stabilize it at some small radius. 

The size of such quantum stabilized vortons is expected to be of the order of the Compton wavelength of the 
fluctuating field - not the astrophysical scales over which the initial string stretched. Nonetheless, vortons are 
expected to be produced during the electroweak phase transition at such high rates that they would naturally evolve 
into acceptable densities of (neutral) dark matter [33]. There are also stringent experimental constraints on the 
density of charged dark matter [44j and thus charged vortons, which may help to rule out certain cosmology models 
that produce such vortons copiously. 

In this section, we investigate a quantum mechanism to stabilize neutral close Nielsen-Olesen type of vortices against 
collapsing to zero radius. The model is that of scalar electrodynamics with the Lagrangian as in eq. ([^, but with 
both the Higgs and the gauge helds now f7(l) valued. The fields of the Nielsen-Olesen string along the infinitely long 
symmetry axis are similar to those in eqs. Q and ([^ with the 2x2 matrices replaced by the identity matrix. They 
are transformed to closed strings by taking the symmetry axis to be of finite length, say 2'kR, and then connecting 
the end points. The symmetry axis thus turns into a circle along which the core of the string is located. When this 
circle is embedded in the x-y plane and centered in the origin, it is appropriate to introduce toroidal coordinates m 
as® 

sinhrcosw sinhrsinw sin cr 

x = R — -, y = R- —r-: z = R- — -.. (32) 

cosh T — cos a cosh r — cos a cosh t — cos cr 

The angle cp G [0, 27r] parameterizes the symmetry transformation along the core circle, i.e. the profiles will not 
depend on this variable. For a constant angle (p, the lines of constant t S [0, oo] are circles that enclose the core of the 
string. Smaller values of r correspond to larger circles that spread out towards infinity on the outside, and approach 
the z-axis on the inside. The inverse 1/t is therefore (in a highly non-linear fashion) related to the distance from the 
string core, similarly to the radius p in the linear Nielsen-Olesen string. Finally, the angle a G [0, 27r] parameterizes 
the points on the fixed r-circles enclosing the string core, i.e. it acts like an azimuthal angle when viewed from the 
core circle. It is therefore natural to give the Higgs field winding by taking a as its phase. Then the profiles are 
predominantly functions of r only, independent of tp, with (probably moderate) modulation in cr, mainly in the vicinity 
of the origin. The origin corresponds to r = 0 and cr = tt, while spatial infinity is approached as r —>■ 0 and cr —> 0. 

Note that the cr-modulation, even if it turns out to be fairly small, cannot be avoided all together: For fixed p and 
T = 0, cr = 0... TT parameterizes the positive z-axis. At each point along this axis the core of the string is seen with 
a different opening angle. For r > 0 but r <C 1, varying cr corresponds to moving along a circle enclosing the string 
core. One side of this circle corresponds to the far distance, while the other side is close to the origin. Obviously, 
different cr values refer to distinct phyiscal regions and the fields must hence depend on this angle. In contrast to the 
linear Nielsen-Olesen string, the profile functions must hence depend on two coordinates, which complicates matters 
significantly. To be explicit, an appropriate ansatz is 

Tl 

® = i'/(o’,'r)e“‘^ and A =—T]{a,T) g{a,T)ea. (33) 


^ See the cited reference m for a visualisation of these coordinates. 
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FIG. 4: The classical energy (left panel) and the energy per unit length of the core circle (right panel) for a closed Nielsen-Olesen 
string as functions of the core radius. The stars indicate the numerically computed cases. Data are presented for g — 
and ms = mw, for which the straight string has energy per unit length 6.28 m^. 


The metric factor rj = coshr — cosct, as well as the inverse coupling constant and radius, have been introduced to 
ensure that / = 1 and 5 = 1 correspond to a pure gauge configuration. Recall that the profiles / and g depend on 
the choice of the dimensionless radius R = mwR- With these preliminaries the classical energy becomes 


E 


dr I dcr sinhr-I ^ ( 9 ^ 5 ) 2 +R ?7 {dr ff + + {d^ ? 7 ^ (l -, (34) 


mw 5" Jo J-n [ Rr] 

subject to the boundary conditions 


w 


/(0,r) =/(27r,r) , 5 ( 0 , r) = 5(27r, r) , and f{a,oo)=g{a,oo) = 0. 


(35) 


Furthermore, Neumann boundary conditions for cr 0, 27r, 


5r/(cr, 0) = drg{<7, 0) = 0 , 


(36) 


ensures that any solution of the field equations with finite total energy will approach / = 5 = 1 as ct —>■ 0 or ct —>■ 27r 
when T = 0 . 

A detailed description of the relaxation method for constructing the actual field configuration numerically con¬ 
strained to a prescribed value of the radius parameter R can be found in ref. |28j . Here we just review the results 
for the minimal energy in figure The left panel essentially reveals an almost linear relation between the minimal 
energy and the radius. More interesting is the energy per length of the torus, which is shown in the right panel. The 
first observation is a saturation at large radii. This is nothing but the energy per unit length of the infinitely long 
string, i.e. the Nielson-Olesen limit. Reproducing this known result is a strong consistency check of our numerical 
treatment. Interestingly, the saturation occurs already for radii only slightly larger than the Compton wave-length 
of the gauge field. This suggests that energy loss from radiating gauge and/or Higgs fluctuations is indeed irrelevant 
for simulations of networks of closed strings that refer to cosmological length scales. For smaller radii, the interaction 
between opposite sides of the torus leads to a considerable drop in the energy per unit length. As a consequence, 
the total energy of the torus configuration decays more than linearly at small core radii i? —>■ 0 , which means that 
it is classically unstable against shrinking to a point. Of course, this instability is expected from Derrick’s theorem, 
according to which higher order derivative terms (as in the Skyrme model[45]) would be needed to stabilize a localized 
object in three space dimensions. 
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Of course, quantum mechanically such an instability could represent a violation of Heisenberg’s uncertainty principle 
according to which field configurations localized within an arbitrarily small region must have large energy and/or 
momentum. To investigate this question, we elevate the radius to quantum variable and treat it as a collective 
coordinate. In a first step we take R = R{t) to be time dependent and then identify its conjugate momentum. Since 


the Cartesian coordinates are fixed in terms of R, this implements a time dependence of the fields via eq. (32). Once 
A is time dependent, its temporal component can no longer vanish but will be induced via Gaufi’ law. As a result, 
it will be proportional to the time derivative R = ■ Thus the Gaufi’ law augments the field configuration by the 

additional profile function q via 


R 

eAo(x,t) = —q{<7,T). 

Putting things together yields the action functional for the unstable mode 


S= dt 


i u{R) R^ - E{R) 


(37) 


(38) 


The energy functional E{R) is that of eq. (34) while the explicit expression for the mass functional u{R) is fairly 
complicated and given in detail in ref. [28j . Here it suffices to mention that u{R) contains linear and quadratic terms 
in the induced profile q. The profile functions / and g are determined from minimizing E{R). In turn they provide 
source terms for q when applying the variational principle to u{R). This procedure is equivalent to solving the Gaufi 
law constraint [28] . The resulting mass functional is shown in the left panel of figure as a function of the radius R. 


The action functional, eq. (38) yields the canonical momentum (operator) p and, via Legernde transformation, the 
(quantum) energy of the unstable mode as 


p = u{R) 


dR 


and 


E = 


+ E{R). 


(39) 


dt 2u{R) 

Figure shows that u{R) is huge for radii larger than the Gompton wave-length of the gauge boson. For a fixed 


momentum, the back-reaction of the classical profiles / and g due to the induced kinetic term in eq. (39) is of the 
order 0{u~^). This effect can thus be safely omitted and these profile functions are reliably obtained by applying the 
variational principle to E{R) alone. As mentioned earlier, q is subsequently obtained from varying u{R), and, since 


initial profiles / and g need not be re-computed, the parameter functions in eq. (39) are then completely determined. 
The corresponding Hamilton operator can then be constructed by imposing canonical commutation relations® for p 
and R. The ground state (g.s.) energy eigenvalue can easily be estimated by recalling that the uncertainty principle 
establishes its lower bound in the form 


This total energy of the torus string is displayed as a function of R in the right panel of figure Obviously a local 
minimum at non-zero torus radius R exists, suggesting the quantum stabilization of the closed Nielsen-Olesen string. 
As for the infinitely long charged string, the length scale at which stability occurs is quite small and it is unclear 
whether it leads to significant cosmological implications. 


® The metric factors in the curvilinear coordinates lead to ordering ambiguities which require particular care. 
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FIG. 5: The mass functional for the unstable mode (left panel) and the lower boundary for its quantum energy (right panel). 
Also shown are the classical energies. 


VIII. CONCLUSION 

This review has mainly focused on computing the fermion contribution to the vacuum polarization energy per unit 
length of an infinitely long straight string in a simplified version of the electroweak standard model. The approach is 
based on the interface formalism, for which the analytical properties of scattering data are essential. Obstacles that 
arise in a naive treatment from the non-trivial structure of the string configuration at spatial infinity are circumvented 
by choosing a particular subset of gauges. Numerically the vacuum polarization turns out to be small and positive 
in the regime in which the one-fermion loop approximation is reliable. Hence, there is no quantum stabilization of 
the (uncharged) string. However, it also turns out that a heavy fermion doublet can stabilize a non-trivial string 
background if such fermions are trapped at the string core to give it a non-zero fixed charge per unit length. The 
resulting configuration is dominated by the Higgs field. Since any additional variational degree of freedom can only 
lower the total energy, the embedding of this configuration in the full standard model, with the C/(l) gauge field 
included, will also yield a bound object, at least if the mixing between this heavy fermion and the known standard 
model fermions can be ignored. Numerically binding sets in at m « 300 GeV, which is still within the range of energy 
scales at which the standard model is expected to provide an effective description of the relevant physics, and also 
within the range to be probed at the LHC. Light fermions would contribute only weakly to the binding of the string, 
since their Yukawa couplings are small. We emphasize that this cosmic string configuration is a novel solution in a 
model that is very closely related to the standard model of particle physics. At a similar scale, though through a 
different mechanism, quantum effects may also prevent neutral closed strings from collapsing. 
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